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Abstract 

C/D 

O 

^^ackground: Investigation of the functioning of the brain in living systems has been a major effort amongst scientists and 

(~l|n'pHiral practitioners. Amongst the various disorder of the brain, epilepsy has drawn the most attention because this 
Oh! 

'-disorder can affect the quality of life of a person. In this paper we have reinvestigated the EEGs for normal and epileptic 
^y^atients using surrogate analysis, probability distribution function and Hurst exponent. 

^^esults: Using random shuffled surrogate analysis, we have obtained some of the nonlinear features that was obtained by 
^;^\ndrzejak et al. [Phys Rev E 2001, 64:061907], for the epileptic patients during seizure. Probability distribution function 
O^hows that the activity of an epileptic brain is nongaussian in nature. Hurst exponent has been shown to be useful to 

o 

(^characterize a normal and an epileptic brain and it shows that the epileptic brain is long term anticorrelated whereas, the 

o 

!riormal brain is more or less stochastic. Among all the techniques, used here, Hurst exponent is found very useful for 
^^Haracterization different cases. 

C^onclusions: In this article, differences in characteristics for normal subjects with eyes open and closed, epileptic subjects 
during seizure and seizure free intervals have been shown mainly using Hurst exponent. The H shows that the brain 
activity of a normal man is uncorrelated in nature whereas, epileptic brain activity shows long range anticorrelation. 
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1 Background 

The brain is a highly complex and vital organ of a human body whose neurons interact with the local as well as the 
remote ones in a very complicated way [TK^IISIHI . These interactions evolve as the spatio-temporal electro magnetic 
field of the brain, and are recorded as Electroencephalogram (EEG) [H[T1[51[B]. Though the detail link between EEGs 
and the underlying physiology is not well understood, the former is widely used for detection and prediction of 
epilepsy, localization of epileptic zone and characterization of the pre and post-ictal [HEllI] using linear and nonlinear 
analysis techniques [Tl[Sl[71[Sl[ni[ini[II] ■ Though mainly nonlinear methods have been applied to predict the onset of 
epileptic seizure and localizing epileptic regions, limited progress has been achieved so far Even some negative 
results have also been reported like linear measures are better than nonlinear measures [121113] , seizure is not a low 
dimensional process [13], it lacks determinism [T51[roi [5]. etc. Hence finding proper analysis techniques is also one of 
the main issues and experts try out different analysis tools for characterizing the normal and diseased brain states, 
especially the epileptic brain. 

In 2001, Ralph G. Andrzejak. et al. and later some other authors [TTiriB] have analyzed five sets of EEG signals [TH] 
each set containing 100 epochs to study the determinism in the brain dynamics for five different physiological and 
pathological conditions. Sets A and B are for normal persons with eyes open and closed respectively and recorded 
extracranially. Sets C and D were recorded intracranially from the hippocampal formation which was 
nonepiletogenic of the opposite hemisphere of the brain and from within the epileptogenic zone of an epileptic 
patient during seizure free intervals respectively. Set E was recorded intracranially from the epileptic zone during 
seizure. The details of the experiments and the conditions have been described in Rcf [1 . R.G. Andrzejak, et al. [T] 
had shown that the normal healthy subject with eyes closed and open shows stochastic behavior using amplitude 
adjusted Fourier transform surrogate analysis where discriminating statistics were the effective correlation dimension 
and nonlinear prediction error whereas, using delay vector variance discriminating statistics, significant nonlinear 
determinism was shown in the same subject |17j . So two conflicting results were obtained for the same subject using 
nonlinear methods. In the case of epileptic patients during seizure and seizure free intervals, determinism was shown 
using two different methods [T[|17j though other studies show lack of determinism for different epileptic patients 
during seizure [ni[ini[I3[IS] . 

On the other hand, characterization of EEGs by scaling properties of the signal is also a major area of research 
interest [8l[9l[Tni[2TJ[22j[23l[24l[25l[26l[27]. Power spectral exponent has been used to characterize the different subjects 
with different physiological conditions [8l[9l[24l[25] and the same exponent has also been used to estimate the 
correlation dimension {Dcorr) 0- Fractal dimension and hurst exponent have also been used to characterize the 
EEGs [26l[27]. Hence a number of experts prefer scaling properties to characterize EEG for different physiological 
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Figure 1: surrogate for data Set A: (a) Variation of Dcorr with m. • for (b) Original data; and A for (c) Surrogate 
data. 

and pathological conditions [5]. 

In this paper, we have reinvestigated the EEG data studied in Refs. [TJ[T71[TH] by random shuffled surrogate analysis 
using Dcorr sls discriminating statistics in order to find determinism in the signal [281I291I5D] and the results have 
been compared with earlier analyses [TJ[T7] . Probability distribution function shows a difference between normal and 
epileptic brain states and this has been discussed in latter Section. Finally, we have quantified the five different 
physiological brain states by Hurst exponent (H) which has been estimated using R/ S analysis |31| . 

2 Results and Discussion 

2.1 surrogate analysis 

Surrogate analysis determines the dynamics in the time series: whether it is governed by stochastic or deterministic 
process [MUMUSO]- 

The surrogate data has been generated by Random Shuffled (RS) surrogate method, in which the signals were 
shuffled randomly so that the probability distribution is same but the temporal correlations are destroyed 
Dcorr which givcs us a measure of the complexity has been estimated for both the original and the surrogate data of 
the data sets A, B, C, D, and E respectively. Fig[IJa) shows that the Dcorr increases with the same trend for both 
the original [Fig[IJb)] and the surrogate data [Fig[ljc)] for the normal persons with eyes open [Set A]. A similar 
trend is observed in the case of the persons with eyes closed, i.e., for Set B [Fig[2I^a)]. This shows that the brain 
activity of a normal person is stochastic in nature agreeing with the analysis by Andrzejak, et al. [T]. For an epileptic 
patient we found a different behavior during seizure free intervals and during seizure activity. The EEG signals, at 
seizure free state, recorded both from the hippocampal formation of the opposite hemisphere of the brain [Fig[2{b) 
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Embedding Dimension (m) 

Figure 2: — • — for original data and — * — for surrogate data for data Set B, C, D and E 

Set C] and within the epileptogenic zone [Fig[2Uc) Set D] show almost same trend in the increase in Dcorr for the 
surrogate data as well as for original data, except for a small separation at higher embedding dimensions, which may 
be due to a high dimensionality of the system. During seizure activity [Fig[2Ud) Set E], Dcorr saturates with 
embedding dimension indicating low dimensional deterministic dynamics and these results agrees well with previous 
analyses [T1[T7]. For sets A-D, as there is no saturation in Dcorr at higher embedding dimension [Fig[lja) and 
Figs[2l^b)-(d)] and hence it is difficult to estimate actual Dcorr- In Ref [1] Dcorr was computed based on quasiscaling 
regions, but such an estimation is very much dependent on the variations of the time- frequency-energy 
characteristics rather than any nonlinear dynamics. Hence this may be inadequate to characterize epilepsy or 
diseased brain states for clinical application [33| . 

2.2 probability distribution functions 

As we have observed from the surrogate analysis that nonlinear dynamics is responsible for epileptic patients during 
seizure, we have compared the probability distribution function (PDF) of a normal case and an epileptic person 
during seizure. The PDF for sets A and E have been shown in Figs[3]Ja) and (b) respectively. Fig[31[a) shows that for 
a normal healthy person with eyes open, the PDF is Gaussian in nature, whereas for epileptic patients during 
seizure, it is nongaussian [Fig[3]^b)] signifying an intermittent nonlinear effect. But for other three cases this feature 
is not so clear. So we feel that the PDF may also be useful to differentiate a brain activity of an epileptic patient 
during seizure from other state. 
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Figure 3: Typical PDF for normal person with eyes closed (up) and for epileptic patient during seizure (bottom). 
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Figure 4: Typical R/S vs lag n plot. 
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Figure 5: Average H with standard deviation error bar. Hurst exponent for hundred time series and average H are 
represented by: □ for normal men with eyes open [set A]; • for normal men with eyes closed [set B]; A and * are 
for the epileptic patients during seizure free interval from two different locations [ set C and D];V for the epileptic 
patients during seizure [set E]. 

2.3 Hurst exponent 

Since one of the major emphasis of epilepsy investigation is to predict their occurrence, it is necessary to know how 
the data is correlated. We have carried out a study of the Hurst exponent (H) which has been estimated using 
Rescaled range analysis (R/S). This method was proposed by Hurst and well established by Mandelbrot, and 
Wallis ini]. For a given set of data series, R/S is defined as [STlfSl] : 

R{n) max{0,Wi,...,Wn) - mm(0, Wi, W„) 



(1) 



Here Wk ^ xi + X2 + + ... + Xk — kX{n), where X, 6*^(71), and n are respectively the mean, variance, and time lag 
of the signal. The expected value of the R/S scales like cn^ as n — > cxo, where H is called the Hurst exponent, and 
can be estimated from the slope of typical plot vs lag (n). For a given signal, we divided the data into 
nonoverlapping blocks of equal length and R/S has been calculated using the Equation [1] and the average value of 
R/S has been plotted as a function of lag in a log — log plot as shown in Fig [4] and estimated H from the slope of the 
curve. For random data H = 0.5, while H > 0.5 for the data with long range correlations, and H < 0.5 indicates the 
presence of long-range anticorrelation or antipersistency in the data. 

The estimated average Hurst exponent (< H >) with an error bar of 100 epochs for all the five EEG data sets (viz. 
A— E) have been shown in FigO The solid box (□) and bullet (•) show the < H > for sets A and B respectively. 
For a normal person with eyes open (set A), the average H (< iJ >) w 0.47 whereas, for the data set B, i.e., for a 
normal person with his eyes closed, < H >« 0.41. 

The < H 0.5 for a normal man with eyes open indicates that the signals are uncorrelated over long time scales 
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Figure 6: Hurst exponent for hundred time series and average H are represented by: □ and — for set A; "and — for 
set B; A and — for set C; * and — for set D;V and for set E. 

signifying stochasticity of the normal brain. But with eyes closed state, decreases inH {< H >w 0.41) may be due to 
the imposition of some extra constraint, which may influence the system towards an antipersistent state. The < H > 
for epileptic patients are shown by up triangle (A) for set C; star (★) for set D; and down triangle (V) for set E 
respectively. < H >« 0.34 and 0.29 for the EEGs recorded at the hippocampal formation and epileptic zone for the 
seizure free intervals of epileptic patients respectively and during seizure (data set E), we get the lowest H 
(< H 0.19). The H for epileptic patient during seizure and seizure free intervals show anticorrelation which may 
be due to epileptiform discharges during seizure free intervals indicating that a large discharge is always followed by 
a small one. The physiology behind the epileptiform discharge is due to the chronic dysfunction or "defect" in the 
epileptic brain, i.e., the epileptic brain is not normal even during seizure free time |35] . Though the hippocampus 
was nonepileptogenic for these subjects, its H is still less than a normal person which may be due to its participation 
in secondary, nonautonomous epileptic processes initiated by the epileptic zone [Tj. The wide dispersion in H for the 
EEG recorded from the epileptic zone in seizure free intervals (*) indicates that the epileptiform discharges are 
intermittent probably due to the chronic presence of abnormal epileptogenic tissues |36[l37j . EEG recording during 
seizure may not be economical and hence it may be better to locate epileptic zone by recording EEG during seizure 
free intervals. These analyses show the possibility of detecting the onset of the seizure state from the time dependent 
Hurst exponent estimated during the transition from normal to seizure state |35j . Fig [5] shows the Hurst exponents 
for all the epochs of five different sets that have been discussed above. 



7 



3 Conclusion 



In this paper we have reinvestigated the EEG data of normal and epileptic subjects to get an insight into the brain 
dynamics at different imposed and diseased conditions using RS surrogate analysis, PDF and H exponents. From 
these analysis we have found that RS and PDF may be useful to find a broad difference between normal and 
epileptic subjects but not helpful for constrained and seizure free intervals. Whereas, using H exponent, we have 
obtained differences in characteristics for normal subjects with eyes open and closed, and epileptic subjects during 
seizure and seizure free interval. The H shows that the brain activity of a normal man is uncorrelatcd in nature 
whereas, epileptic brains show long range anticorrelation. 
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